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A  GRAPHICAL  SIMILARITY  MEASURE  FOR  TIME  SERIES  MODELS 

Benjamin  Kedem 

University  of  Maryland,  College  Park 

Abstract 

The  oscillation  in  time  series  can  be  represented  by  counts  of  axis- 
crossings  in  the  series  and  its  differences.  These  counts  are  called 
higher  order  crossings  and  display  a  monotone  property  whose  rate  of  in¬ 
crease  discriminates  between  processes.  Only  very  few  of  these  counts 
are  needed  for  effective  discrimination  as  shown  by  plots  of  higher  order 
crossings  obtained  from  real  and  simulated  data. 
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1.  Introduction 


The  purpose  of  this  note  is  to  introduce  a  graphical  device  useful 
as  a  measure  of  similarity  or  as  a  goodness  of  fit  criterion  for  hypothe¬ 
sized  time  series  models.  It  is  based  on  the  actual  oscillation  observed 
in  time  series  as  depicted  by  axis-crossings  and  higher  order  crossings. 
Higher  order  crossings  (HOC)  are  axis-crossings  of  differenced  time  series 
and  are  closely  linked  to  the  spectral  content  of  the  series.  In  fact 
under  the  Gaussian  assumption,  to  which  we  shall  adhere,  HOC  determine 
the  finite  dimensional  distributions  up  to  a  scale  parameter  given  that 
the  mean  is  zero.  The  main  advantage  of  HOC  is  that  they  are  easily  ob¬ 
tained  from  an  observed  series  and  that  only  very ^ few  of  them  are  needed, 

{'■ 

as  the  discriminatory  power  in  HOC  usually  diminishes  with  their  order. 

Higher  order  crossings  in  time  series  discrimination  were  discussed 
in  Kedem  and  Slud  (1981),  (1982), where  a  certain  goodness  of  fit  criterion 
is  suggested.  Here  however  the  emphasis  is  on  a  graphical  device  rather 
than  a  single  test  statistic.  This  graphical  method  may  be  shown  useful 
in  answering  the  question  "Does  a  given  time  series  oscillate  as  a  certain 
hypothesized  model?*  Some  examples  with  real  and  simulated  data  demon¬ 
strate  the  use  and  potential  of  this  method. 


2.  Definition  and  Properties  of  HOC 


Let  { } ,  t  =  0,±l,...,  be  a  zero  mean  stationary  Gaussian  process 
with  correlation  function  and  spectral  distribution  function  F.  The 
higher  order  crossings  of  order  k,  is  the  number  of  axis-crossings 


yk  *z  yk-1z 


where  V  is  the  difference  operator 


Vzt  *  zt  '  zt-i 
k  k-1  0 

and  V  Z  =  V(V  Z  ).  Since  V  is  the  identity  operator,  D  is  the 

t  L  i  N 

actual  number  of  axis-crossings  by 

A  more  precise  definition  is  constructed  in  terms  of 
1  Vk-17  >  0 

=  c  ,  k  =  1,2, . . . 

0,  otherwise 


_  Xtk>  * 

0,  otherwise. 


We  define  formally  the  higher  order  crossings  of  order  k  by 


n  -  +  +  a  k 

k,n  2  N 


In  practice  we  only  deal  with  finite  series  and  some  care  must  be 
taken  in  regard  to  the  differencing  of  finitely  long  records.  Each  time 
a  finite  series  is  differenced  an  observation  is  lost  so  that  if  k  higher 
crossings  are  desired  the  time  index  t  =  i  is  given  to  the  k'th  or  a 


later  observation  and  we  shall  follow  this  rule. 


Consider  the  following  example.  Suppose  the  series  7  is  given  by 
-1  -4  -3  .6  -.3  -.8  1.2  -1.3  -.4  -1.1. 


In  order  to  derive  the  first  three  HOC  D.  ,,  D„  , ,  D~  ,  we  have 

i  t)  ^  >  D  J  j  6 


t 

1 

2 

3 

4 

5 

6 

z 

-1 

.4 

-  .3  .6 

-  .3 

-.8 

1.2 

-1.3 

-.4 

-1.1 

x(1) 

0 

1 

0  1 

0 

0 

1 

0 

0 

0 

vz 

1.4 

-  .7  .9 

-  .9 

-.5 

2.0 

-2.5 

.9 

-  .7 

x(2) 

1 

0  1 

0 

0 

1 

0 

1 

0 

_2 

7  Z 

-2.1  1.6 

-1.8 

.4 

2.5 

-4.5 

3.4 

-1.6 

(3) 

X 

0  1 

0 

1 

1 

0 

1 

0 

and 

Dl,6=2’ 

°2, 6 

=  4’  D3,6  =  4 

and  we 

say  that  D^, 

°2*  °3 

were 

derived 

from  a  series  of  length  6.  It  is  seen  that  counting  symbol  changes  in  the 
clipped  binary  series  is  equivalent  to  counting  axis-crossings  in  the  cor¬ 
responding  series.  Thus  ^  = 2  since  there  are  2  symbol  changes  in 

or  equivalently  2  axis-crossings  in  Z^,...,Z^.  It  is  in¬ 
teresting  to  observe  that  while  D  is  the  number  of  crossings,  D„  ,  is 

o  2.  i  o 

the  number  of  peaks  and  troughs  in  Z,,...,Z, . 

1  D 

The  fact  that  the  correlation  structure  in  { }  is  completely  deter¬ 
mined  by  the  sequence 


“VJ-i 


is  due  to  the  basic  relation 


Pj_  =  cosi 


(1) 


whose  extension  is  given  in  (2)  by  considering  the  first  correlation  in  V  Z 


U +pi((>U}  + 


[k)-2pi[H]  +  -  +  ^2pk 


From  (1;  and  (2)  we  can  determine  p^,  P2>  •••  recursively. 

Kedem  and  Slud  (1981).)  (2)  can  also  be  written  as 


cos(ui)  (sin  IjO))  dF(w) 


/rED,  ,  .  I  c°s(“)<s 

I _ k+1,  N\  _  _J[ _ _ _ 

In-1/"  f 

(sin 

—TV 


which  is  referred  to  as  the  higher  order  crossings  spectral  representa¬ 
tion  and  relates  the  sequence  (ED^  to  F.  This  spectral  representation 
has  been  recently  studied  in  Kedem  (1984)  in  some  detail,  where  it  was 
shown  that 


ED1,N  5  ED2,N 


<  (N-l)  . 


r  ,  oo 

It  follows  that  (ED.  1 .  .  is  a  monotone  increasing  and  bounded  sequence 

J  .N  j  =  l - — - - 

and  therefore  converges.  It  can  be  shown  that 


<A  j 


where  w*  is  the  highest  point  in  the  support  of  F,  that  is,  the  highest 
frequency  In  the  spectrum.  Therefore 


EDl,N  5  ED2,N 


I"  (N-l) 


ic 

In  practice  most  series  contain  noise  and  u )  =  tr.  When  tu  =  it  the  in¬ 
equalities  in  (4),  (6)  are  strict.  It  was  shown  in  Kedem  and  Slud  (1982) 
that  when  n  is  in  fact  included  in  the  support  of  K  then  the  actual 


.  1  it.  W  .1  ■  *  •  i  *  V  .  j  -  '  »  E 


.  «  .  •  k  0  "j 


5 


sequence  (D.  }  is  monotone  increasing  provided  N  is  sufficiently  large. 

Now  it  is  this  monotone  property  of  HOC  which  we  would  like  to 

utilize.  More  precisely,  we  wish  to  demonstrate  that  the  initial  rate 

of  increase  in  the  D.  as  j  ■><x>  is  rather  fast  and  carries  discrimina- 

J » “ 

tory  information  as  different  processes  display  different  rates  of  in¬ 
crease  in  the  D.  for  small  j  and  fixed  N.  As  i  increases  this 
discriminatory  potency  decreases  rapidly  and  similar  rates  are  displayed 

by  many  different  series.  In  this  sense  only  very  few  I).  are  needed 

J  >  ” 

in  achieving  effective  discrimination. 

It  should  be  noted  that  from  (2),  (3),  the  information  contained 

GO 

in  {ED^  j^}j_i  is  equivalent  to  that  in  {p  }  and  consequently  that  in  F, 
the  normalized  spectral  distribution  function.  But  the  information 


manifested  by  the  ED.  ^  presents  a  fferent  angle  of  view  and  another 
way  for  the  interpretation  of  stationary  time  series.  HOC  provide  an¬ 
other  dimension  in  time  series  analysis,  which  addresses  the  oscillatory 
content  in  time  series  directly  by  counts  of  visual  features. 


3.  An  Approximation  to  the  Variance  of  HOC 


As  we  shall  restrict  our  attention  to  the  rate  of  increase  in  HOC 
it  is  necessary  to  determine  beforehand  the  variance  of  D  ,  k=l,2,...  . 

K ,  In 

In  general  this  variance  is  a  function  of  4'th  order  orthant  probabili¬ 
ties  wniah  are  not  easily  accessible  (see  Reed  (1983)).  It  has  been  ob- 
served  however  that  in  many  cases  the  dependence  in  the  {d^  }  for  low  k 

is  rather  weak  as  expressed  by  very  low  correlations,  unless  the  original 

series  {z^}  displays  extremely  large  absolute  correlations.  Thus  as  a 

(k) 

first  approximation  we  suggest  treating  {d^  ' }  as  a  binary  Markov  chain. 
This  simplification  leads  immediately  to  a  computable  expression  for 
Var(Dk  N)  which,  as  we  shall  see,  provides  a  rather  close  approximation 
to  the  true  variance. 

Define 


p(x[k)  =  1 


,(k) 

Vj 


1) 


Pj(k)  =  Corr(VkZt  ,  V1^.). 


Then  from  (1)  or  (2)  we  have 


X) 

J 


(k)  =  1  +  1  _,_-l,(k-l) 


„  sin  p  . 
TT  J 


(7) 


(k) 

Observe  that  because  of  (1),  is  precisely  given  by  (2). 

(k) 

parameters  which  specify  (d^  s  are  thus  given  by 


The  two 


P<k>  =  P(d^k)=l)  =  i-x{k) 


(8) 


1  ?.\(k)  +  \(k) 

(k)  _  ,  ,(k)  _  ,  ,(k)  J  2 

v  =  P(d  -  1  |  d  -  1) -  (k) 

t  t  )  2(1  -  X  ;) 


(9) 


7 


Therefore  using  standard  arguments  from  the  theory  of  Markov  chains  (see 

Karlin  and  Taylor  (1975),  Ch.  2)  we  have  the  approximation  owing  to  the 

fact  that  D.  x.  =  d^  +  ...  +  df^ , 
k ,  h  2  N 


/ex  /ex  o„(k)  (k)  k  )  .  j 

Uar'V»>  *  W-1*!'  1  +  'TS"‘>'Vk'N 


where 


\,N  =  <« 


J  /  ..(k)  -  (k)  \  N-ll  /  .  . 

[ 1  -  (-7^-)  \/a-'’  1 


(k)  (k) 

and  q  =1  -  p  .  Usually  V  is  negligible  compared  with  N  and  may 

K  >  IN 

be  omitted.  Thus,  a  useful  approximation  to  Var(D  )  is 

K  j  N 

—  -1 

,,  .  ...  (k)  (k)  1  -  2£^k)  +  v(k)  ,. 

Var^Dk  "  (N-I)p  9  - - (kT~  ' 

From  a  given  sequence  p  ,p  ,...,p  we  compute  Var(D  ),  k=l,.. 

1  2.  K.  9  IN 

K-l  as  follows. 

fkl 

(ij  Obtain  "  ^  from  (2),  k  = 0, 1 , . . . , K-l . 

-n+pi((2kM2k)}-P  4(2k)+(2k)}+-+‘-*>v- 

(k)  _  \ -  r  lLMi_  __k_-21-L  k-3/  J 

P1  '  ''  '  '  /  2k\  /  2k \  k 

( +"’ +  (_1)  2Pk 


(10 

Then  is  given  by 


„<k>  -  -  t  +  2p«  -  2p<k+1)(J-„<k>). 


(ii)  Obtain  A^k\  A.^  from  (7),  k-l,..., K-l. 

,00  1.1.  -1  (k-l)  ,(k)  _  ] 

A1  =  2  +  if  Sin  P1  *  A2  -  2  ^ 


,  (k)  _  1  1  .  -1  (k-l) 

A2  -  2  +  n  sln  1 1 2 


8 


(iii)  Obtain  P(k),  v(k),  from  (8),  (9),  k  =  l,...,K-l 


P(k)  =  1-  A<k), 


i_9A(k)+A(k) 
(k)  =  1  2A1  +  A2 

2(1  -Ajk)) 


(iv)  Substitute  in  (10)  to  obtain  Var(D  ),  k  =  l,...,N-l. 

K,  IN 

The  above  algorithm  results  are  compared  in  what  follows  with  actual 
estimates  in  the  estimation  of  Var(Dj  iqqq) *  •  =  1,...,6  for  several 

series.  To  apply  the  algorithm  one  needs  to  supply  as  input  only  p  .,..., p^ 
and  N.  (A  computer  program  which  takes  as  input  p^,...,p_,N  and  gives  as 
output  E(D.  M)  and  Var(D.  )  is  available  upon  request.)  When  an  hypothe- 
sized  model  is  entertained  these  parameters  are  usually  available  as  is 
the  case  for  example  when  the  hypothesized  process  is  white  noise  or  any 
specific  autoregressive  moving  average.  The  experimental  results  were  ob¬ 
tained  from  100  independent  realizations  of  size  N  = 1000  each  and  are 
summarized  in  Tables  1,  2,  3.  Note  that  ED^  j^qqq  available  exactly 
and  does  not  depend  on  the  Markov  assumption.  The  tables  also  compare 

ED.  , with  their  estimates  as  a  check  on  the  simulation  results.  It 
J , 1000 

is  seen  that  ED.  . _ and  ED.  are  very  close  throughout.  In  pariicu- 

j,1000  j,1000 

u 

lar,  {Var  D.  2  obtained  from  (10)  agrees  well  with  actual  estimates 

obtained  from  these  100  realizations.  We  note  that  throughout  our  simu¬ 
lation  V  in  (10)  ranged  in  value  from  0.663  to  1 . 2 A 7  but  most  often 
k ,  N 

its  value  was  around  0.9,  which  is  very  small  compared  with  N=  1000.  in 
specifying  a  model  we  use  the  notation  as  in  Box  and  Jenkins  (1970).  for 
example,  AK(1)  with  parameter  (p  means  7.  -  ^  +  uc,  but  MA(I)  witli 

parameter  (J  means  7.^  =  u  —  Ou^  j,  where  u^  are  N(0,0^)  independently  dis¬ 
tributed.  It  should  be  noted  that  aLl  our  results  concerning  HOC  are  scale 


free  so  that  the  actual  magnitude  of  is  not  important 
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FIRST  ORDER  MODELS 


Series 

j 

EDj,iooo 

“j.iooo  ,Var  “j.iooo1'5 

from  (10) 

(^r  aj,1000 
from  100 
realization: 

White 

1 

500 

497 

15.81 

15.96 

Noise 

2 

666 

666 

13.15 

13.63 

3 

7  32 

732 

12. 16 

12.53 

4 

769 

770 

11.57 

11.49 

5 

794 

795 

11.18 

11.05 

6 

813 

814 

10.82 

10.00 

AR(1) 

1 

333 

333 

15.  74 

15.58 

2 

580 

582 

14.00 

15.03 

<f>  =0.5 

3 

685 

686 

12.72 

13.11 

4 

740 

740 

12.00 

12.08 

5 

774 

774 

11.49 

11.23 

6 

797 

796 

11.09 

10.00 

MA  (1) 

1 

369 

371 

14.17 

14.59 

2 

553 

551 

12.44 

12.27 

8  =  -0.  5 

3 

640 

639 

11.91 

12.48 

4 

694 

693 

11.66 

12.22 

5 

732 

731 

11.44 

11.40 

6 

759 

759 

11.22 

11.02 

Table  1. 

Comparison  of 

fVar  Dj,1000J 

from  (10)  and 

the  estimate 

obtained  from 

100  independent 

realizations 

each  of  size 

1000  given  in  the  last  column.  ED.  , ___  are  rounded  to 

j,1000 

the  nearest  integer. 


SECOND  ORDER  MODELS 


Series 

j 

KD 

j , 1000 

EDj , 1000 

{Var  ,1000^ 

r  “j.lOOO1 

from  (10) 

f  rom  1 00 
realizations 

AR(2) 

1 

333 

333 

12.00 

10.83 

<j>^  =  0.  75 

2 

3 

460 

562 

460 

562 

11.87 

12.40 

11.07 

12,58 

4u  =-0.5 

4 

646 

649 

12.65 

13.58 

£ 

5 

709 

12.25 

12.57 

6 

751 

751 

11.79 

11.91 

MA(2) 

1 

342 

343 

18.33 

17.76 

0  =-0.5 

0 

3 

648 

782 

649 

784 

17.66 

13.82 

18.06 

13.71 

0„  =-0.8 

4 

835 

836 

5 

857 

858 

10,24 

9.56 

6 

868 

869 

9.  74 

8.71 

AR(2) 

1 

424 

425 

9.64 

9.67 

=  0.4 

2 

3 

484 

536 

485 

537 

9.38 

10.29 

9.13 

10.81 

=-0.7 

4 

594 

594 

11.27 

12.72 

L 

5 

651 

652 

11.87 

12.02 

6 

702 

12.04 

11.34 

Table  2. 


Comparison  of  {Var  D.  ^qqq) 2  from  (10)  with  the  estimate 

obtained  from  100  independent  realizations  of  size  1000 

given  in  the  last  column,  ED.  , are  rounded  to  the 
°  ] , 1 000 

nearest  integer. 
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A.  Probability  Limits  for  HOC 

The  asymptotic  distribution  of  as  can  be  found  under 

standard  moment  conditions  without  recourse  to  the  Markov  assumption 
which  we  only  needed  in  getting  an  approximation  to  the  variance  of  HOC. 
For  any  finite  moving  average  series  (z^l,  {d^ }  is  m-dependent  for  some 
finite  m  which  implies  (Diananda  (1953))  that  the  asymptotic  distribu¬ 
tion  of  D  as  N  is  normal.  More  generally,  since  (Z  }  is  Gaussian 

K  j  IN  C 

the  condition  J  |p^|  <  00  implies  the  asymptotic  normality  of  1)^  ^  for 
any  finite  k.  The  proof  of  this  fact  uses  a  method  suggested  by  Malevich 
(1969)  and  Cuzick  (1976)  and  was  used  in  Kedem  (1980),  Ch.  7,  in  deriving 
the  asymptotic  distribution  of  D1  .  The  extension  to  D  is  immediate 
and  will  not  be  reproduced  here.  It  follows  that  when  p.  is  absolutely 

K 

summable  D  is  asymptotically  normal  with  approximate  variance  given 

K  >  IN 

by  (10),  and  hence  approximate  95%  probability  limits  for  0j_  ^  are 

(N-l)p(k)  ±  1 . 96{Var  DR  ^  (12) 

(k) 

where  Var(D  )  is  given  by  (10),  and  p  is  as  in  (8).  (12)  provides 

a  graphical  means  for  assessing  the  similarity  or  dissimilarity  between 
processes  as  expressed  by  their  oscillatory  information  by  way  of  higher 
order  crossings.  Several  examples  follow. 

Table  A  gives  the  probability  limits  (12)  for  three  processes  whose 
data  are  given  in  an  appendix  in  Priestley  (1981).  It  is  seen  that  the 
HOC  obtained  from  these  records  with  N  =  A50  fall  well  within  the  proba¬ 
bility  limits.  The  graphical  display  of  the  HOC  and  the  limits  (12)  ol 
these  data  are  given  in  Figure  1,  and  show  graphically  the  similarity  of 
the  observed  data  and  the  hypothesized  processes.  At  the  same  time  the 
figure  shows  the  dissimilarity  or  the  distance  of  the  particular  AK(2) 


and  ARMA(2,2)  under 

consideration 

from  white 

noise,  and  also  the  dis- 

similarity  of  these 

AR(2)  and  ARMA(2,2).  In  particular,  the  rates  of 

increase  in  the  D. 

J» 

displayed  by  the  three 

processes  are  different 

and  serve  as  a  fast 

discrimination 

feature. 

Rounded 

95% 

Rounded 

probability 

limits 

EUj,450 

U  j  ,450 

203 

246 

225 

221 

282 

317 

299 

299 

313 

345 

329 

328 

330 

361 

346 

348 

342 

372 

357 

355 

351 

380 

365 

363 

a)  White  noise 

:  Z  =  u 

Table  4.  Observed  in  records  given  in  an  Appendix  in 

PriestLey  (1981).  The  ARMA(2,2)  series  c)  contains  a 
missing  observation  which  was  replaced  by  0.0. 

~  N(0, 1),  independent. 
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Figure  1. 


c)  Observed  D.  from  an  ARMA(2,2)  series  and 

J  ,  A  jU 

their  theoretical  limits  (12). 


HOC  representation  for  the  three  series  in  Table  4.  In  each 

case  the  theoretical  limits  (12)  capture  the  correct  curve 

obtained  from  observed  D.  ,rn  from  the  indicated  series. 

3,450 

(The  dark  lines  mark  the  probability  limits.) 


Our  next  example  concerns  a  HOC  diagnostic  chock  for  time  series 
models  as  applied  to  the  temperature  series  C  in  Box  and  Jenkins  (1976) 
Our  first  step  is  to  see  whether  the  centered  series  and  its  first  dif¬ 
ference  are  close  to  being  white  noise.  A  quick  glance  at  Figure  2-(n) 
shows  that  the  series  and  its  difference  are  far  from  being  white  as 
their  HOC  do  not  fall  within  the  white  noise  limits  with  N=212.  On 
the  other  hand,  the  HOC  of  the  fitted  model  (V%t  =  0.82VZt  ^  +  u  ) 
suggested  by  Box  and  Jenkins  (1971),  p.  293,  arc  well  within  the*  proba¬ 
bility  limits  (12)  as  seen  from  Figure  2-(b).  The  residuals  HOC  are 
given  in  Figure  2-(c)  and  are  shown  to  fall  well  within  the  white  noise 
limits  except  for  ^j_2  which  still  very  close  to  being  "in".  Thus 
based  on  the  oscillatory  properties  of  the  differenced  series  and  the 
corresponding  residuals  it  is  seen  that  the  fitted  model  is  reasonable. 


j 

957  limits  for 
white  noise 

95%  limits  for 
AR(1) ,  0  =  0.82 

_ 

D.  from 

J  1 212 

Z  -  Z 

D j ,212  from 

vzt 

D.  from  the 

J  ,212 

residuals 

-  0.82VZtl 

1 

91 

118 

28 

54 

5 

48 

89 

2 

129 

153 

98 

125 

48 

105 

139 

3 

143 

166 

129 

153 

105 

137 

155 

4 

152 

173 

144 

166 

137 

153 

168 

5 

158 

178 

152 

173 

153 

166 

172 

6 

162 

182 

158 

178 

166 

172 

171 

Table  5.  Observed  HOC  from  the  temperature  series  C  in  Box  &  Jenkins 

(1976)  and  probability  limits  for  two  models.  (N  =  212).  The 
HOC  from  VZ^  fall  within  the  AR(1)  limits.  The  residuals  HOC 
suggest  that  the  residuals  very  much  resemble  white  noise. 
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a)  HOC  from  Z  -Z  and  VZ 
t  t 

fall  outside  the  white 


noise  limits. 


b)  HOC  from  VZ  fall  inside  the 
limits  (12)  under  the  model 
VZt  =  0.82VZ  j  +  u  . 
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c)  HOC  of  the  residuals  series  resemble  very 
much  HOC  of  white  noise. 


Figure  2.  HOC  analysis  of  the  tempe-nture  series  C  in  box  &  Jenkins 
(1976)  with  N  =  212.  (The  dark  lines  mark  the  probability 
limits. ) 
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Our  last  example  concerns  periodic  data  plus  noise-  In  many  appli¬ 
cations  a  major  problem  is  whether  an  observed  series  is  made  of  noise 
only  or  of  signal  plus  noise.  The  oscillation  in  periodic  signals  pius 
noise  is  usually  different  from  the  oscillation  displayed  by  white  noise 
and  cen  be  captured  very  well  by  HOC.  Consider  the  sunspot  series  given 
and  analyzed  in  Anderson  (1971).  It  is  well  established  that  this  series 
contains  several  significant  periodic  components,  and  from  our  point  of 
view  it  is  interesting  to  measure  the  deviation  of  this  series  from  pure 
white  noise  as  depicted  by  higher  order  crossings.  By  appealing  to  the 
limits  (12)  with  N  =  155  we  can  see  from  Figure  3  that  the  series  does 
indeed  contain  a  "signal"  and  the  hypothesis  of  white  noise  is  rejected. 
That  is,  the  sunspot  series  does  not  oscillate  as  white  noise. 


i 

95%  probability  | 

limits  for  white  noise 

D  .  ,  c  c  f  ram 
.1*155  _ 

Zt-Z 

l 

65 

89 

31 

2 

93 

113 

36 

3 

103 

122 

81 

4 

110 

128 

107 

5 

114 

131 

1 1  9 

6 

117 

134 

121 

Table  6.  Observed  HOC  from  the  sunspot  series 

(N *  155)  in  Anderson  (1971)  and  proba¬ 
bility  limits  for  white  noise.  The 
hypothesis  (  f  noise  is  rejected. 


5.  Some  Power  Calculation 

In  many  cases  and  D ^  alone  are  capable  of  discrimination  between 

time  series  models.  However,  there  are  many  examples  in  which  the  later 

D  do  the  job  while  D^,  fail.  In  order  to  illustrate  this  point  we 

consider  several  examples  of  autoregressive  moving  average  models  where 

the  hypothesis  is  that  of  white  noise.  At  the  same  time  the  examples 

indicate  the  power  of  our  graphical  method  obtained  from  50  independent 

series  each  of  length  N  = 450.  In  each  case  we  bring  in  the  results  of 

10  realizations  in  the  form  of  binary  row  vectors  of  length  6.  bach 

realization  produces  D,  /r_,...,D.  ,  and  when  a  certain  0.  , falls 
r  1,450  6,450  J.450 

outside  the  probability  limits  (12)  for  the  white  noise  case,  it  is  in¬ 
dicated  by  "1"  for  success.  Otherwise  this  is  indicated  by  "0".  Thus 
110011  means  that  D, ,  D„,  Dc ,  D  fall  outside  the  limits  (12)  for  the 
white  noise  case,  while  D^,  fall  inside  these  limits.  Clearly,  for 
series  which  are  not  white  noise  we  expect  to  see  many  l's.  The  closer 
the  process  is  to  white  noise  the  harder  is  the  discrimination  problem 
and  the  ratio  of  0's  increases.  It  is  seen  from  Figure  4  that  there 
are  series  which  are  far  from  being  white  but  their  actual  axis- 
crossings  behave  as  those  of  white  noise.  On  the  other  hand,  the 
difference  in  oscillation  is  captured  by  higher  order  crossings  as  seen 
in  particular  from  (g)  and  (i).  The  power  is  seen  to  increase  with  more 
significant  parameters  as  is  well  expected. 
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a )  White  noise  b)  AR(1),  <1>  =  0  -  05  c)  MA(1 ),  0-0.1  d)  AR(1),  <p-.2 
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power  =  0.4  power  =  0.9 
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0  110  0  0 

0  0  110  0 
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-0.15  g)  ARMA(2,2),  = -.2 

0,  =  .2,  0 2  =  .1 ,  power  =  0.88 
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h)  ARMA(1,1),  <J>,  =  .1,  0,  =-.l  i) 
power  =  0.86 
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Figure  4.  Ten  typical  binary  vectors  obtained  from  IK  =  in 

relation  to  the  white  noise  limits  from  (12).  The  power  was  esti¬ 
mated  from  50  such  binary  vectors.  1  means  corresponding  I). 
is  on  I  s  i  do  l  he  limits:.  Otherwise  0  is  recorded. 
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